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Time series of satellite data provide unparalleled information on the response of vegetation to climate variability. 
Detecting subtle changes in vegetation over time requires consistent satellite-based measurements. Here, the 
impact of sensor degradation on trend detection was evaluated using Collection 5 data from the Moderate 
Resolution Imaging Spectroradiometer (MODIS) sensors on the Terra and Aqua platforms. For Terra MODIS, 
the impact of blue band (Band 3, 470 nm) degradation on simulated surface reflectance was most pronounced 
at near-nadir view angles, leading to a 0.001 -0.004 yr -1 decline in Normalized Difference Vegetation Index 
(NDVI) under a range of simulated aerosol conditions and surface types. Observed trends in MODIS NDVI over 
North America were consistent with simulated results, with nearly a threefold difference in negative NDVI trends 
derived from Terra (17.4%) and Aqua (6.7%) MODIS sensors during 2002-2010. Planned adjustments to Terra 
MODIS calibration for Collection 6 data reprocessing will largely eliminate this negative bias in detection of 
NDVI trends. 

© 201 1 Elsevier Inc. All rights reserved. 


1. Introduction 

Satellite data capture changes in the land surface, atmosphere, and 
ocean on diurnal to decadal time scales. Large and remote areas can be 
monitored using satellite data, and at nearly 40 years, the lengthening 
time series of satellite observations offer insight into climate impacts on 
human and natural systems. In particular, satellite-based studies of 
changes in vegetation phenology from seasonal, annual, and decadal- 
scale temperature variability have received considerable attention in 
recent years (e.g., Goetz et al., 2005; Myneni et al., 1997; Wang et al., 
2011). Detecting subtle changes in vegetation characteristics requires 
long time series of data from well-calibrated sensors (Slayback et al., 
2003) and careful attention to uncertainties associated with atmospheric 
correction (Nagol et al., 2009) and sensor degradation. 

Data from the Moderate Resolution Imaging Spectroradiometer 
(MODIS) sensors on board the NASA Terra (2000-present) and Aqua 
(2002-present) satellites provide a unique opportunity to assess the 
consistency of satellite-based measurements. Data products from both 
sensors are generated using consistent algorithms, and MODIS standard 
products support a range of science applications (e.g., Esaias et al., 1998; 
Justice et al., 2002; King et al., 2003). Recently, several studies have 
reported inconsistent results from the Terra and Aqua MODIS sensors, 
including estimates of ocean chlorophyll-a concentrations (Djavidnia 
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et al., 2010) and aerosols (Levy et al., 2010). Such discrepancies may im- 
pact the interpretation of time series or even the ability to detect tem- 
poral trends, especially for MODIS products that utilize data from both 
Terra and Aqua MODIS sensors. 

Here, the impact of sensor degradation on trends in MODIS normal- 
ized difference vegetation index (NDVI) was analyzed over boreal forest 
and tundra cover types in North America, biomes previously identified 
as undergoing rapid climate-driven changes in productivity (de Jong 
et al., 2011; Pouliot et al., 2009). Degradation coefficients for each 
MODIS spectral band were calculated, and the time-varying trends in 
MODIS degradation were propagated through a simplified version of 
the MODIS atmospheric correction process. The simulated NDVI trends 
were then compared to observed trends in NDVI from the Terra and 
Aqua MODIS sensors over boreal North America. This study indicates 
that negative NDVI trends of up to — 0.004 yr - 1 can be attributed solely 
to degradation of the visible and near-infrared (NIR) bands aboard the 
Terra MODIS instrument. Therefore, careful consideration of sensor 
degradation may be warranted when interpreting decadal scale trends 
in land, ocean, and atmospheric constituents based on Collection 5 
(C5) Terra MODIS data. 

2. Data and methods 

2.1. MODIS top-of-atmosphere reflectance 

Sensor degradation is a common occurrence as satellite instruments 
age in the harsh space environment. In theory, sensor calibration using 
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onboard apparatus can correct for sensor degradation, but in practice, 
the solar diffuser, stimulation lamps, and stability monitors needed for 
onboard calibration also change during mission life. A useful secondary 
source of information is to monitor reflectance trends with “pseudo-in- 
variant” desert targets, which appear to be radiometrically and spectral- 
ly stable during the satellite era (Cosnefroy et al., 1996). 

MODIS C5 top-of-atmosphere (TOA) reflectance data were collected 
over stable desert sites (Wu et al., 2008), endorsed by the Committee on 
Earth Observation Satellites (CEOS) (http://calval.cr.usgs.gov/sites_ 
catalog_map.php), to estimate degradation of the Terra (2000-2010) 
and Aqua (2002-2010) MODIS sensors over time. Level lb TOA data 
for MODIS Bands 1 (red), 2 (NIR), and 3 (blue) were extracted from 
16-day repeatable orbits over the site to maintain consistent sensor 
viewing angles. TOA reflectances were normalized using site-specific 
bi-directional reflectance function (BRDF) coefficients determined 
from the first three years of Terra and Aqua observations (Wu et al., 
2008). The two sides (mirror sides 1 and 2) of MODIS scan lines had 
different performance and were treated separately. The BRDF normal- 
ized TOA data time series for each view angle (frame), spectral band, 
and mirror side were fitted using a second-order polynomial to 
estimate degradation for days since launch (DSL): 

Degradation = a 0 + a A ■ DSL + a 2 ' DSL 2 ( 1 ) 

where, a 0 , ai and a 2 are regression coefficients. To estimate the impact 
of planned changes in sensor calibration for Collection 6 (C6) MODIS 
reprocessing (scheduled for late 2011), the C5 degradation coefficients 
were adjusted using the ratio of C5:C6 TOA reflectance for the Libya 4 
site. Prototype C6 TOA data were generated over selected targets to 
evaluate the impact of calibration changes on MODIS data products. 
For the first time, C6 will incorporate calibration information from 
pseudo-invariant targets. A preview of the potential for C6 data to cor- 
rect the bias in Terra MODIS products is provided here based on C6 cal- 
ibration approach II, presented in May 201 1 at MODIS Science Meeting 
(http://modis.gsfc.nasa.gov/sci_team/meetings/201105/). C6 data for 
the entire MODIS time series will not be available until at least 2012. 


Therefore, the potential for continued degradation in C5 Terra data 
was also estimated by extending the fitted degradation curves to 2013. 

2.2. Simulation of atmospheric radiative transfer and atmospheric correction 

The time-dependent degradation coefficients for C5 Terra MODIS 
data were used to simulate the impact of sensor degradation on trends 
in NDVI over boreal North America during 2000-2010. The simulation 
included three major steps. First, 1 1 -year time series of TOA reflectance 
were generated for forest and tundra cover types under a range of aero- 
sol conditions, adjusting TOA values for sensor degradation over time. 
The proposed simulation approach was similar to the methods used to 
generate the MODIS surface reflectance product (MOD09, (Vermote et 
al., 2002; Vermote & Kotchenova, 2008)), with several simplifying as- 
sumptions. The 6S vector atmospheric radiative transfer model (Version 
1.1, (Kotchenova et al., 2006)) was used. To improve computational 
efficiency, look-up tables (LUTs) of spherical albedo (p), path reflec- 
tance (r 0 ), and transmittance (7) were constructed. These LUTs were 
based on 15 aerosol optical depth (AOD) values (0.01-2.0), three aero- 
sol models (continental, background desert, and biomass burning), and 
small angular intervals (5° steps for solar and view zenith angles and 
10° steps for relative azimuth angle). The LUTs were used to relate 
TOA reflectance (r) and surface reflectance (r s ) (Liang, 2004): 


Mid-summer surface reflectance inputs (r s ) for different cover types 
were derived from Aqua MODIS NDVI data (MYD13C2, (Huete et al., 
2002)). NDVI values were binned in increments of 0.1, and mean spec- 
tral reflectance values were estimated for each NDVI bin. Using Eq. (2), 
TOA reflectance (r) for each NDVI bin was then estimated under differ- 
ent aerosol conditions using the LUTs described above. Finally, TOA re- 
flectance values (r) for each year were multiplied by the degradation 
coefficients estimated for C5 Terra MODIS data to simulate the TOA 
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Fig. 1. Degradation coefficients of red (B1 ), NIR (B2), and blue (B3) top-of-atmosphere reflectance for Terra and Aqua MODIS C5, and Terra MODIS future C6 data. Degradation co- 
efficients were calculated separately for data from Mirror Side 1 (MSI) and Mirror Side 2 (MS2). Additional VZAs are shown in Fig. SI. 
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Fig. 2. Projected degradation of Terra MODIS C5 red (Bl), near-infrared (B2), and blue 
(B3) top-of-atmosphere reflectance in years 2011-2013 (dashed lines) for four 
near-nadir VZAs. 


reflectance (r ) observed by the sensor for each view zenith angle (VZA) 
and spectral band. 

In the second step, simulated Terra MODIS surface reflectance values 
(r s ) were derived from the degraded TOA reflectance values (r). A 
simplified Dense Dark Vegetation algorithm (Kaufman et al., 1997) 
was used to retrieve AOD. Two simplifications were made: 1) direct 
use of Band 3 surface reflectance, rather than estimating Band 3 reflec- 
tance using other MODIS bands (Kaufman et al., 1997), since C5 Terra 
MODIS degradation impacts were minimal for longer wavelengths, 
and 2) a single estimate of AOD at 550 nm from Band 3, thereby ignor- 
ing the spectrum dependence of AOD, because the aerosol model was 
known in the simulations. Given AOD, Eq. (2) can be solved to obtain 
surface reflectance (r s ): 

f = f ~ r ° (3) 

s 7 + (r-r 0 )p 


Finally, time series of NDVI values were calculated using the estimat- 
ed red and NIR surface reflectance values (r s ) (Huete et al., 2002) for 
each NDVI bin, year, and set of atmospheric conditions (AOD and 
aerosol model). 

Two additional sets of simulations were developed to characterize 
the sensitivity of the results to specific input parameters for AOD and 
surface reflectance. In Scenario 1, the trend in the MODIS Terra AOD 
product (MOD08_D3) over land (Levy et al., 2010) was used instead 
of simulated AOD trends based on Band 3 degradation for the atmo- 
spheric correction procedure. This scenario tested the sensitivity of 
the results to estimates of blue band degradation, since the MODIS 
AOD product over land uses the blue reflectance in combination with 
other spectral information to estimate AOD (Remer et al., 2005). Sce- 
nario 2 examined the impact of sensor degradation in the red and NIR 
bands on simulated trends in NDVI. In this second scenario, the impact 
of blue band degradation was simulated while holding inputs for red 
and NIR reflectance constant. This scenario examined the degree to 
which degradation of the Terra MODIS Bands 1 and 2 partially compen- 
sates for Band 3 degradation over time. 

2.3. Monte Carlo simulation 

While the radiative transfer modeling estimated the effect of sensor 
degradation on individual NDVI time series, it did not address the issue 
of how such degradation might confound the extraction of vegetation 
trends from a large (e.g. continental) set of data. Of particular concern 
was estimating how the degradation-induced errors influence the re- 
sults of detection of NDVI trends, given the stated uncertainty of the 
NDVI product. To explore this issue, a set of Monte Carlo simulations 
were used to examine the sensitivity of NDVI time series trend detec- 
tion to artifacts from sensor degradation. Time series of C5 NDVI from 
peak growing season conditions in 2000-2010 were used as inputs to 
the Monte Carlo simulation. Each time series had an “actual” linear 
NDVI trend (T 0 ) of varying magnitude ( — 0.02 yr _1 to +0.02 yr -1 ), 
plus a linear degradation. Three degradation rates (no degradation 
and degradation of 0.001 and 0.003 NDVI yr - 1 ) were used. Indepen- 
dent and identically distributed Gaussian errors were then added to 
the NDVI values. The standard deviation of the errors was set 0.015, 
based on MODIS validation studies over AERONET sun photometer 
sites (http://modis-sr.ltdri.org/). The sensor degradation was then 
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Fig. 3. Simulated Terra MODIS C5 NDVI (solid lines) and AOD (dashed lines) time series for near-nadir observations (VZA= — 11) of tundra (left) and forest (right) during 
2000-2010 with low (0.05, top), moderate (0.2, middle), and high (0.5, bottom) AOD. Simulated NDVI and AOD time series for all VZAs are shown in Fig. S2. 
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Table 1 

Annual trends in simulated Terra MODIS C5 NDVI during 2000-2010 over tundra (NDVI = 0.65) based on moderate aerosol loading (AOD = 0.2). Results were calculated separately 
for each VZA and mirror side (MS) using three different aerosol models. Negative VZAs refer to pixels from the beginning of a scan line to nadir; positive VZAs indicate pixels from 
nadir to the end of a scan line. Stars “*” indicate no statistically significant trend. Values shown here are trends multiplied by 10 3 . 
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applied to generate the final NDVI time series. The time series were fit 
with a linear trend, using a t - test and 95% confidence interval to identify 
statistically significant linear trends. 

Three variables (bias, ratio of false negatives, and ratio of false 
positives) were calculated for each set of Monte Carlo simulations. The 
bias in trend detection was estimated using the mean of all statistically 
significant trends, p{T Detected ), and the actual trend, T 0 : 

Bias = ^Detected) — To (4) 

M) 


series were fit with a linear trend, using a t - test and 95% confidence 
interval to identify statistically significant linear trends. For pixels 
with significant NDVI trends in more than one compositing period, 
only pixels in which all significant trends had the same slope direction 
(positive or negative) were retained for the final map. A map of forest 
disturbances between 1990 and 2000 was used to exclude areas of 
fire damage or timber harvest predating the MODIS record (Masek 
et al., 2008), and likely disturbance areas during the MODIS era were 
excluded using a threshold of NDVI change between sequential obser- 
vations (>0.2). 


The ratio of false negatives was estimated as the fraction of all time 
series (N) in which detected trends did not have the correct positive or 
negative signs (n incorre ct) or the trends did not pass the significance test 

( Hno_sig ) • 

r ^ incorrect "T ^Inosig /£ -\ 

F n = jv 

The ratio of false positives was estimated as the fraction of all statis- 
tically significant trends ( n sig ) in which statistically significant trends 
had the opposite signs ( n inC orrect ): 


2.4. MODIS NDVI trends for North America 

For comparison with the simulated impacts of sensor degradation 
described above, trends in C5 Terra and Aqua MODIS NDVI for boreal 
North America during 2002-2010 were mapped. High quality NDVI 
values from the C5 MODIS 500 m vegetation indices product 
( MOD 1 3 Al /MYD 1 3 Al , (Huete et al., 2002)) were analyzed. Data were 
selected for four 16-day composites during the growing season 
(starting on day of year 177, 193, 209, and 225 for Terra MODIS; day 
of year 185, 201, 217, and 233 for Aqua MODIS) to minimize impacts 
of vegetation phenology, and time series analyses were conducted 
separately for each compositing period. Only pixels with high quality 
NDVI observations in at least 6 out of the 9-year record were included 
in the analysis based on NDVI pixel reliability estimates. NDVI time 


3. Results 

3.1. Degradation of top-of-atmosphere reflectance 

Sensor degradation varied by spectral band, view angle, and mirror 
side (Fig. 1, SI). Overall, Aqua MODIS Bands 1-3 showed <2% change 
in TOA reflectance between 2002 and 2010, while Terra MODIS data 
exhibited substantial degradation, particularly at shorter wavelengths 
(Fig. 1, SI ). Terra MODIS Band 3 had the most pronounced degradation, 
with a decrease of nearly 7% at near-nadir view angles over the past 
decade. By 2013, degradation of C5 Terra MODIS Band 3 may exceed 
10% since launch (Fig. 2). Band 3 is not used directly for calculating 
NDVI; instead, degradation of Band 3 TOA reflectance over time will im- 
pact the calculation of surface reflectance in other spectral bands 
(Kotchenova et al., 2008; Vermote et al., 2002) and higher-level 
MODIS products, including NDVI. Planned calibration changes for C6 
Terra MODIS data reprocessing will largely eliminate trends in TOA 
reflectance observed in the C5 time series (Fig. 1). Changes in TOA 
reflectance for prototype C6 Terra MODIS for 2000-2010 were <2% for 
most cases. 

3.2. Surface reflectance simulations 

Degradation of the Terra MODIS blue band resulted in decreasing 
NDVI trends (2000-2010) for simulated tundra and forest cover types 
under a range of aerosol conditions (Fig. 3). Changes in NDVI arose 
from the decreasing trend in retrieved AOD from degraded blue TOA re- 
flectance. Three aerosol models were tested, yet the choice of aerosol 
model had little impact on retrieved AOD values in the simulations. 


Table 2 

Annual trends in simulated Terra MODIS C5 NDVI for Scenario 1 based on observed Terra MODIS AOD decline over land ( — 0.003 yr - 1 ). Results were calculated separately for each 
VZA and mirror side (MS) using three different aerosol models. Negative VZAs refer to pixels from the beginning of a scan line to nadir, whereas the positive VZAs refer to pixels 
from nadir to the end of a scan line. Values shown here are trends multiplied by 10 3 . 


View zenith angle 

-51 

-42 

-32 

-28 

-18 

-11 

-3 

9 

17 

31 

43 

52 

MSI 

Continental 

-1.2 

-1.0 

-0.8 

-0.8 

-0.5 

-0.5 

-0.5 

-1.3 

-0.4 

-0.6 

-0.3 

-0.5 


Desert 

-1.3 

-1.1 

-0.8 

-0.8 

-0.6 

-0.6 

-0.5 

-1.4 

-0.5 

-0.6 

-0.4 

-0.6 


Bio burning 

-1.2 

-1.0 

-0.7 

-0.7 

-0.5 

-0.5 

-0.5 

-1.3 

-0.4 

-0.5 

-0.3 

-0.5 

MS2 

Continental 

-1.2 

-0.8 

-0.8 

-1.1 

-0.6 

-0.5 

-0.7 

-0.7 

-0.6 

-0.5 

-0.5 

-0.5 


Desert 

-1.4 

-0.9 

-0.8 

-1.1 

-0.7 

-0.6 

-0.8 

-0.8 

-0.7 

-0.6 

-0.6 

-0.7 


Bio burning 

-1.2 

-0.8 

-0.7 

-1.0 

-0.6 

-0.6 

-0.7 

-0.7 

-0.6 

-0.5 

-0.5 

-0.6 


D. Wang et al. / Remote Sensing of Environment 119 (2012) 55-61 


59 


Table 3 

Annual trends in simulated Terra MODIS C5 NDVI for Scenario 2 in which blue band degradation was simulated while red and NIR bands were kept constant. Results were calculated 
separately for each VZA and mirror side (MS) using three different aerosol models. Negative VZAs refer to pixels from the beginning of a scan line to nadir, whereas the positive 
VZAs refer to pixels from nadir to the end of a scan line. Stars “*” indicate no statistically significant trend. Values shown here are trends multiplied by 10 3 . 
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For a moderate (AOD = 0.2) and high (AOD = 0.5) aerosol loading, 
retrieved AOD declined by approximately 0.1 between 2000 and 2010. 
Under the low aerosol case, retrieved AOD reached the minimum simu- 
lation value (0.01 ) as the blue reflectance decreased over time. 

Negative trends in NDVI varied by VZA and aerosol loading 
(Table 1, SI ). For moderate AOD, NDVI over forest and tundra declined 
by 0.003 yr - 1 at near-nadir VZA but <0.001 yr - 1 for VZA >30° from 
nadir. Simulations with high AOD had even stronger NDVI declines 
from sensor degradation (>0.004 yr -1 , Table SI). In general, rates of 
declining NDVI were similar for tundra and forest types when compar- 
ing simulations with the same AOD and VZA (Fig. 3, S2). The impacts of 
sensor degradation on time-varying trends in NDVI were further 
elaborated in two additional experiments. In Scenario 1, which was 
based on trends in Terra AOD, declines in AOD (0.003 yr -1 ) and NDVI 
(0.001 yr -1 ) were approximately l/3rd of the simulated declines 
based on degradation of Bands 1-3 (Table 2). In Scenario 2, excluding 
degradation of Terra MODIS red and NIR bands had little influence on 
the simulated NDVI decline at near-nadir VZAs (Table 3) since the 
degradation coefficients of red and NIR bands were similar (see Fig. 1). 

Implementing calibration changes planned for C6 MODIS reproces- 
sing largely eliminated trends in Terra MODIS C5 NDVI from sensor 
degradation over 2000-2010. For moderate AOD over tundra, a limited 
number of off-nadir view angles retained small positive NDVI trends 
(0.001 yr -1 ) in simulations using C6 calibration information (Table 4). 

3.3. NDVI trend detection 

The systematic NDVI decrease from Terra MODIS sensor degradation 
could bias trend detection using C5 data. Based on the Monte Carlo 
simulations, the bias in trend detection could be > 500% when the actual 
NDVI change is small (<0.002 yr -1 ) and the degradation rate is large 
(0.003 yr -1 ). In this case, nearly all detected trends originated from 
sensor degradation rather than actual vegetation dynamics (Fig. 4). 
When degradation of NDVI was simulated at 0.003 yr -1 , almost half 
of simulated NDVI trends of + 0.006 yr - 1 could not be detected or 
were detected with the opposite sign. Even a degradation-based NDVI 
decline of 0.001 yr -1 resulted in >50% false negative detections of all 
NDVI trends between —0.002 and + 0.004 yr -1 . Overall, the Monte 
Carlo simulations demonstrate that the combined impact of sensor 


degradation and product uncertainties have a large impact on correct 
identification of trends in simulated C5 Terra MODIS NDVI. 

3.4. Terra and Aqua NDVI trends 2002-2010 

Both Terra and Aqua MODIS indicated regions with statistically 
significant negative (browning) and positive (greening) NDVI trends 
during 2002-2010 (Fig. 5). More than 95% of trends detected using 
2002-2010 C5 Terra and Aqua data for North America exceeded the 
estimated threshold for trend detection (0.004, Section 3.3). However, 
there was a threefold difference in the area of negative NDVI trends 
from Terra compared to Aqua, with large contiguous regions where 
Terra data indicated small negative NDVI trends. NDVI trends over bore- 
al forest and tundra vegetation were small, with the mode of Aqua NDVI 
trend of + 0.008 year - \ Monte Carlo simulations (Section 3.2) suggest 
that these subtle changes fall within the range of time series trends that 
may be impacted by the combination of sensor degradation and product 
uncertainties. The combination of subtle changes in NDVI and compara- 
ble impacts from sensor degradation resulted in substantial differences 
in the spatial patterns of NDVI trends from Terra and Aqua. For pixels 
with statistically significant NDVI trends from both sensors, Terra 
NDVI trends were systematically smaller than Aqua NDVI trends during 
2002-2010 (Fig. 6). The mean difference between Terra and Aqua NDVI 
trends was — 0.002 yr -1 (Fig. 6). This value falls within the range of 
simulated declines in Terra MODIS NDVI from sensor degradation. 

4. Discussion and conclusions 

Degradation of the Terra MODIS blue band resulted in negative NDVI 
trends over tundra and boreal forest vegetation, with greater impacts at 
the near-nadir VZAs favored in the MODIS reflectance and NDVI prod- 
ucts. Simulated impacts of sensor degradation were similar to observed 
trends in C5 Terra MODIS data for North America. Previous studies 
based on Advanced Very High Resolution Radiometer (AVHRR) and 
Landsat data have identified greening and browning trends of 
0.005-0.01 NDVI yr - 1 in the boreal and tundra regions of North Amer- 
ica (de Jong et al., 201 1 ; Pouliot et al., 2009). Monte Carlo simulations in 
this study suggest that up to half of all actual NDVI trends of this magni- 
tude could be missed by Terra MODIS when both sensor degradation 


Table 4 

Annual trends in simulated Terra MODIS C6 NDVI data during 2000-2010 over tundra (NDVI = 0.65), based on anticipated calibration changes in C6 reprocessing. Simulations were 
calculated for moderate aerosol loading (AOD = 0.2) at each VZA and mirror side (MS) for three different aerosol models. Negative VZAs refer to pixels from the beginning of a scan 
line to nadir, whereas the positive VZAs refer to pixels from nadir to the end of a scan line. Stars “*” indicate no statistically significant trend. Values shown here are trends mul- 
tiplied by 10 3 . 


View zenith angle 

-51 

-42 

-32 

-28 

-18 

-11 

-3 

9 

17 

31 

43 

52 

MSI 

Continental 

* 

* 

1.3 

1.2 

* 

* 

* 

* 

1.1 

* 

* 

* 


Desert 

* 

* 

1.3 

1.2 

* 

* 

* 

* 

1.1 

* 

* 

* 


Bio burning 

* 

* 

1.1 

1.1 

* 

* 

* 

* 

1.0 

* 

* 

* 

MS2 

Continental 

* 

* 

1.6 

* 

* 

* 

* 

* 

1.0 

* 

* 

* 


Desert 

* 

* 

1.6 

* 

* 

* 

* 

* 

1.0 

* 

* 

* 


Bio burning 

* 

* 

1.4 

* 

* 

* 

* 

* 

0.8 

* 

* 

* 


60 


D. Wang et al. / Remote Sensing of Environment 119 (2012) 55-61 



Actual NDVI trend 


Fig. 4. Monte Carlo simulation results for the combined impacts of NDVI degradation and 
NDVI uncertainty (0.015) on trend detection. Detection bias (top), percentage of missing 
detections (middle), and percentage of false detections (bottom) were summarized for all 
statistically significant linear trends for different actual NDVI trend levels. 


and product uncertainties are considered. Initial simulations using data 
processed with prototype C6 algorithms indicate no significant trends 
associated with sensor degradation for near-nadir VZAs. Thus, it is 
expected that C6 will largely remedy the issues described in this 
paper. Although MODIS data over the North American boreal forest 
and tundra were used in this study, the impacts of sensor degradation 
are not limited to this region because sensor degradation impacts all 
Terra/MODIS C5 TOA data and subsequent processing for land and at- 
mosphere products. In addition, degradation of Terra MODIS C5 data 
impacts ongoing efforts to inter-calibrate data from current and previ- 
ous generations of Earth observing satellites to generate long time se- 
ries of NDVI (e.g., Fensholt & Sandholt, 2005; van Leeuwen et al., 
2006). Based on temporal trends in C5 Terra MODIS data identified in 
this study, we recommend that users rely on Aqua MODIS data for 
time series analysis of vegetation changes and inter-calibration of 
MODIS, AVHRR, and the next generation of polar orbiting satellites 
(e.g., the NPOESS Preparatory Project— NPP) until C6 data are available 
for the Terra MODIS time series. 


Studies of subtle trends in vegetation cover may be particularly 
sensitive to artifacts from Terra MODIS sensor degradation. Many land 
products select near-nadir observations where sensor degradation 
was most pronounced, and the small negative bias in Terra MODIS 
NDVI from sensor degradation may be more easily confused with real 
browning trends over vegetation than cloud contamination or other 
sensor artifacts. Given the non-linear declines in Band 3 TOA reflec- 
tance, studies that compare annual or seasonal anomalies in C5 Terra 
MODIS NDVI may also amplify artifacts from sensor degradation. 
Many other applications will be less sensitive to artifacts from sensor 
degradation, such as studies that use longer wavelengths or applica- 
tions based on large and persistent changes in surface reflectance asso- 
ciated with land cover change. 

Terra MODIS Band 3 degradation may also account for observed 
discrepancies in AOD from Terra and Aqua MODIS sensors. Levy et al. 
(2010) noted that Terra's AOD record showed a negative trend 
compared to ground-based AERONET sunphotometers while compari- 
sons between AERONET and Aqua AOD showed no trend. Over land, 
MODIS aerosol retrieval is based on three bands, including Band 3 
(Remer et al., 2005). Sensitivity tests suggest that >3% degradation of 
Band 3 is sufficient to explain the trend in Terra MODIS versus AERONET 
comparison. As with studies of vegetation, C5 Aqua MODIS data should 
be used for studies of inter-annual variability in aerosols until Band 3 
degradation can be addressed in C6 reprocessing. 

Sensor calibration is critical for time series analyses. Long time se- 
ries desired by the science community necessitate consistent calibra- 
tion throughout the entire mission lifetime. In the case of MODIS, 
calibration efforts now utilize onboard equipment, monthly lunar ob- 
servations, and pseudo-invariant ground targets to determine time- 
dependent and wavelength-dependent changes in the reflective 
solar bands. Instrument characteristics and on-orbit malfunctions 
also complicate sensor calibration. The observed degradation on 
Terra MODIS is thought to be associated with on-orbit changes in 
the bi-directional reflectance characteristics of the solar diffuser 
from those obtained in pre-launch tests, which is used along with 
measurements by the solar diffuser stability monitor to derive on- 
orbit degradation rates. Following an anomaly with the Terra 
MODIS solar diffuser in 2003, the solar diffuser door was fixed in 
the “open” position to keep the instrument screen permanently in 
place (Xiong & Barnes, 2006). Since that change, the rate of Terra 
MODIS degradation has increased. These issues highlight the impor- 
tance of having multiple lines of evidence for assessing sensor calibra- 
tion, including the use of both on-board mechanisms and invariant 
ground targets. 

Terra MODIS C5 data remain an important source of information 
for a range of science applications. However, time series analysis of 



Fig. 5. Annual NDVI trends calculated from 2002 to 2010 Terra/and Aqua/MODIS C5 NDVI products over North America. Terra/MODIS C5 NDVI data indicated largely negative 
(browning) trends, while statistically significant trends in NDVI from Aqua/MODIS were primarily positive (greening). Recent disturbances from fire and wood harvest appear 
white, while water, cropland, and barren land cover types appear black. Inset panels summarize statistically significant NDVI trends over land from each dataset. 
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- 0.01 - 0.008 - 0.006 - 0.004 - 0.002 0 0.002 0.004 0.006 0.008 0.01 

Fig. 6. Difference in 2002-2010 annual NDVI trends (as shown in Fig. 5), for pixels where both Terra and Aqua MODIS indicated statistically significant trends based on C5 data. 
Differences were calculated as Terra-Aqua. Differences were largely negative, indicating that Aqua trends were more positive at each pixel than Terra trends (histogram, inset). 
Overall, pixels detected as significant trends in both datasets represented 34.3% of trend pixels detected by Terra and 44.0% of pixels detected by Aqua. 


vegetation or aerosol properties based directly or indirectly on C5 
Terra MODIS Band 3 data may be sensitive to impacts of sensor deg- 
radation. Errors from sensor degradation remain within the uncer- 
tainties of the retrieval algorithm and instrument characteristics 
(Vermote & Kotchenova, 2008). However, even small trends in sensor 
performance may influence time series analysis. C6 data reprocessing 
will likely alleviate the problems identified in this study, yet MODIS 
data users should consider the potential for sensor degradation to im- 
pact their analysis, given that the aging MODIS sensors will remain 
vulnerable to similar problems in the future. 

Appendix A. Supplementary data 

Supplementary data to this article can be found online at doi:10. 
1016/j.rse.201 1.12.001. 
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